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Long-term (a few days) variation of magnetic helicity injection was calculated 

for 28 solar active regions which produced 47 CMEs to find its relationships with 

C/3 ' the CME occurrence and speed using SOHO /MDI line-of-sight magnetograms. As 

a result, we found that the 47 CMEs can be categorized into two different groups 

by two characteristic evolution patterns of helicity injection in their source active 

regions which appeared for ~0.5-4.5 days before their occurrence: (1) a monoton- 
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ically increasing pattern with one sign of helicity (Group A; 30 CMEs in 23 active 
regions) and (2) a pattern of significant helicity injection followed by its sign rever- 
sal (Group B; 17 CMEs in 5 active regions). We also found that CME speed has 
a correlation with average helicity injection rate with linear correlation coefficients 
Q\ ■ of 0.85 and 0.63 for Group A and Group B, respectively. In addition, these two 

CME groups show different characteristics as follows: (1) the average CME speed 
of Group B (1330kms _1 ) is much faster than that of Group A (870kms _1 ), (2) the 
CMEs in Group A tend to be single events, whereas those in Group B mainly consist 
of successive events, and (3) flares related to the CMEs in Group B are relatively 
more energetic and impulsive than those in Group A. Our findings therefore suggest 
that the two CME groups have different pre-CME conditions in their source active 

H ■ regions and different CME characteristics. 
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Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: pho- 
tosphere — Sun: magnetic topology — Sun: surface magnetism — Sun: evolution 



Introduction 



Coronal Mass Ejections (CMEs) are transient ejections into interplanetary space of as 
much as about 10^ 3 -10 17 g of plasma and embedded magnetic fields from the solar corona 
flGopalswamyl 120061 ) . There have b een many studie s to un derstand CME initiation with sev- 
eral numerical simulation models. IChen fc Shibatal ( 120001 ) showed in their simulation that a 
CME can be triggered by the localized reconnection between a pre-existing coronal field and 
a reconnection-favored emerging magnetic flux system with some observational supports (e.g., 



Fevnman &: Martini Il995l; iJing et a 



van Ballegooiien &: Martens 



1989 



Linker et al. 



Martini Il994t iGaizauskas et al. 



1997 



2004: IWang et all 12004 ). A flux cancellation model (e.g., 



20011). i n agreement with some observations 



Martini Il998l ). suggested that flux cancellation at 



(e.g., 

the neutral line of a sheared arc ade can lead to the for mation, destabilization, and eruption of 
a filament followed by a CME. lAntiochos et al.l ( 119991 ) proposed a model (so-called breakout 
model) in which the reconnection of a sheared arcade with an overlying background magnetic 
field gradually removes a constraint over the sheared arcade so that a CME can occur. A kink in 



stability (e.g- lriood fc Priestlll98~ll:lFan fc GibsoJ2004l JTorok fc Kliemll2005l:lcho et alJboogh 



and a torus instability (e.g., lBatemanlll978l : iTitov fc Demoulinlll999l : iKliem fc Torokll2006l ) by 
an e mergence of a t wisted flux tube were also considered to explain the initiation of CMEs, 



e.g., iFan fc Gibson! (120071 ) performed isothermal MHD simulations of the three-dimensional 
evolution of coronal magnetic fields as a twisted magnetic flux tube emerges gradually into a 
pre-existing coronal arcade. In addition, there are som e other trigger mechanisms: a shear (o r 
twist) motion of footpoints of magnetic arcade s (e.g. 



filament mass drainage f e.g..lLow 



2001 



Mikic et al. 



slow decay of a background magnetic f i eld (e.g.. Ilsenberg et al.lll993l ). a buoyancy force due to 



19881 : iKusano et al.ll2004 . 



Zhou et al.ll2006l ) , and Moreton & EIT waves generated 



by a remote CME (e.g., iBallesterl 120061 ) . However, a precondition and a trigger mechanism for 
CME initiation has not been understood clearly yet. More observational studies are needed to 
determine and/or suggest a most convincing model which can explain the detailed processes of 
CME occurrence. 

After CMEs are initiated at the corona, they depart from the Sun at speeds ranging from 
~20 to more than 3000 km s" 1 (average sp eed of 480 kms" 1 ) m easured from the SOHO /LASCO 
white-light images during solar cycle 23 (jGopalswamyl 120061 ) . To understand how CME speed 
is associated with magnetic properties in CME-productive active regions, there have been some 
studies to find a relationship between CME speed and several magnetic parameters derived from 
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photospheric magnetic fields. iQiu fc Yurchyshynl (]2005h studied 13 CME events a nd found that 
there is a strong correlation between CME speed and total reconnection flux (see Forbes fc Lin 
2000l . and references therein) estimated from p hotospheric magnetic fields with a linear correla- 
tion coefficient (CC) of 0.89. iGuo et all (120071 ) examined properties of photospheric line-of-sight 
magnetic fields in 55 active regions before the onset times of 86 CMEs originating from the 
investigated active regions. They measured four magnetic parameters in the active regions 
(i.e., tilt angle, total flux, total length of strong-field strong-gradient neutral lines, and effective 
distance) and found a moderate correlation between CME speed and effective distance with a 
linear CC of -0.4. 

Magnetic helicity has been recognized as a quantitative and useful measure for global 
complexity and non-potentiality of a given magnetic field system in terms of twists, kinks, and 



inter- linkages of magnetic field lines fjBerger fc Field! Il984l ; iPevtsovl I2008J ; iDemoulin Sz Pariat 



20091 ). However, there have been only a few studies of magnetic helicity that is related to 
solar eruptive phenomena such as filament erupt ions and CMEs . By s tudying a filament in the 
active region NOAA 8375 on 1998 November 4, iRomano et al.l (120031 ) found that the filament 
has a helically twisted configuration (shown by a TRACE 171 A image) and there is a steady 
magnetic helicity injection of 2xl0 41 Mx 2 from the photosphere of the active region into the 
corona over a time span of ~28 hours around the times of four eruptive events involved in 
the filament. From the study of 34 pairs of CMEs and magnetic clouds (MC), ISung et al. 



( 120091 ) showed that the square of CME speed is positively cor r elated with magnetic helicity per 



unit length in the corresponding MC at 1 AU. ISmyrli et al.l ( 120101 ) investigated the temporal 
variation of photospheric magnetic helicity injection in 10 active regions which produced 12 
halo CMEs. They found that there is no unique trend of a short-term (several hours) variation 
of magnetic helicity injection before and after CME occurrence even though a r apid and large 



chang e of magnetic helicity injection was observed in some cases. Recently, IRomano et al. 



( 1201 ll ) examined a filament eruption in the active region NOAA 9682 and reported that the 
filament eruption may be caused by the interaction of two magnetic field systems with opposite 
signs of magnetic helicity. 

In this study, long-term variation (a few days) of magnetic helicity injection is investigated 
for 28 active regions which produced 47 CMEs. The goals of this study are: (1) to find a 
characteristic variation pattern of helicity injection in relation to a pre-CME condition and a 
CME energy build-up process and (2) to carry out a correlation study between CME speed 
and helicity injection rate in the active regions. The rest of this paper is organized as follows: 
data selection and analysis are explained in Section [2J In Section [31 we first describe in detail 
the observational finding of two CME groups classified by two characteristic evolution patterns 
of helicity injection and their different CME characteristics (§ 13. ip . then the relationship of 
CME speed with helicity injection rate for the two CME groups is presented (§ 13. 2p . Finally, 
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we discuss in Section H] that the two helicity patterns can be interpreted by some numerical 
simulation models for CME initiation. 



2. Data Selection and Analysis 

For a statistical study of magnetic helicity injection in active regions related to CMEs, we 
selected 28 active regions which produced 47 CMEs (refer to Table HJfor the detailed information 
of the CME even ts). Out of total 47 CMEs, 42 CMEs were adopted from the CME list of 



Guo et all ( 120071 ) in which the CMEs are identified with their originating active regions. We 
carefully checked the identification of the CME source regions by investigating not only the 
CME position angles with respect to their source regions but also CME-related phenomena in 
the source regions such as soft X-ray flares, EIT brightenings/dimmings, erupting filaments, or 
post-flare arcades. The same method for the CME source region identification was applied to 
add five more CME samples in this study. 

By helicity, we refer to the relative magnetic helicity in the rest of this paper, i.e., the 
amount of helicity subtracted by that of the corresponding potential field state. The injection 
rate of helicity H through the photospheric surface S of the CME-produc tive active regions at 



a specific time t is calculated from the method developed by I Chad ( 120011 ): 



H{t) ~ / G A (x,t)dS, (1) 



s 

where the integrand Ga{x, t) = — 2[v^cT(x,t) ■ A p (x,t)]B n (x,t) is a good proxy of helicity 
flux density (i.e., helicity injection per unit area per unit time) at a point x on S. GA{x,t) is 
calculated using MDI magnetogram data with a spatial resolution of 2" per pixel as follows: (1) 
B n (x,t), the magnetic field perpendicular to S, is derived from Bi(x,t), the MDI line-of-sight 
magnetic field, based on the assumption that the transverse component of the magnetic field 
on S is negligible compared to the longitudinal component (i.e., B n (x,t) = Bi(x, t)/cosi/j(x) 
where ip(x) is the heliocentric angle of the point x). However, to reduce the uncertainty of 
determining B n (x,t), we calculate Ga(x, t) only if x is located within 0.6 of a solar radius 
from the apparent disk center, (2) i>lct(#, t), the velocity of the apparent horizontal motion 
of field lines on S, is determined by the technique of local correlation tracking (LCT) using 
two consecutive 96-min MDI magnetogram data, and (3) A p (x,t) is a specific vector potential 
satisfying: 

h-V x A p (x,t) = B n {x,t), (2) 

V ■ A p (x, t) = 0, (3) 

A p (x,t)-h=0. (4) 
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A p (x,t) is calculated by using the fast Fourier transform (FFT) method. 

Figure []Ja and [Tb present ^lct and A p of AR 10720, respectively, which are plotted as 
arrows on the grayscale image of B n at 14:23 UT on 2005 January 16. The grayscale Ga map 
is also shown in Figure [lb: positive value (i.e., helicity flux density of right-handed sign) of Ga 
is displayed as white tone, while ne gative value (i.e . , helic ity flux density of left-handed sign) of 
Ga is displayed as black tones. See IChae fc Jeongl (120051 1 for the details of the Ga calculation. 



Many previous stu dies have ta 
activ e regions (e.g., 



ten acc o unt of Ga for calcu 



Vloon et al. 



20051 : IJeong fc Chae 



20071 : 



2002a, 



LaBonte et al. 



Nindos et al. 



2007 



2003] 



Lim et al 



ating helicity injection in solar 



Chae et al. 



20071 : 



2004 



Park et a 



Romano et al. 



2008 



201 



3) . Not e 



that although a better proxy of helicity flux density Gg was proposed by IPariat et al.l (120051) . 
howev er, by the comparison between these two methods (Ga a nd Gg) of H calcu lation, I Chae 



(J20071 ) found that their discrepancy is typically less than 10 %. iPark et al.l ( 120101 ) found that 
the relative standard deviation of H calculated from the same method as used in this study 
is around 5 %. We therefore think that the uncertainty of Ga does not significantly affect the 
helicity calculation and conclusion of this study. 

After H is determined as a function of time, the accumulated amount of helicity injection 
AH in the CME-productive active regions is determined by integrating H with respect to time: 



AH(t) 



Hit') dt', 



(5) 



to 



where to is the start time of each MDI magnetogram data set of the active region under 
investigation, to is set as the time when the active region appears or rotates to a position 
within 0.6 of a solar radius from the apparent disk center. As a reference parameter, we also 
derive the total unsigned magnetic flux $ at S: 



$(t) 



\B n (x,t)\dS. 



(6) 



Results 



3.1. Two Characteristic Patterns of Helicity Injection before CMEs 

We have investigated helicity injection through the photospheric surface of 28 solar active 
regions over a span of several days around the times of 47 CMEs which originated from the 
active regions. As a result, it was found that all the CMEs under investigation are always 
preceded by a significant helicity injection of 10 42 -10 44 Mx 2 over a long period (~0.5-4.5 days) 
into the active-region corona through the photosphere. Note that there have been previous 
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reports of a large amount of helicity injection in active regions for ~0. 5 -4.5 days before the 



occur rence of major flares (e.g. , lLaBonte et al.ll2007t iPark et al. 



2008 



20101 ). a filament eruption 



(e.g.. iRomano et al.l 120031 ). and halo CMEs (e.g.. ISmyrli et al.l 120101 ). Furthermore, we found 
that there are two characteristic patterns in the long-term (a few days) variation of helicity 
injection in the CME-productive active regions which appeared for a few days before the CME 
occurrence: (1) a monotonically increasing pattern with one sign of helicity (Group A; 30 
CMEs in 23 active regions) and (2) a pattern of significant helicity injection followed by its sign 
reversal (Group B; 17 CMEs in 5 active regions). 

Figures [21 [3J and H] present the time profiles of AH calculated for the 23 active regions in 
Group A. As shown in the time profiles, each of the 30 CMEs in Group A was preceded by a 
large injection of one sign of helicity. More specifically, the absolute value of AH significantly 
increased at a nearly constant rate, (0.3-110) x 10 40 Mx 2 hr" 1 , over about 0.5-4.5 days and then 
26 out of the 30 CMEs occurred during the monotonically increasing phase. The other 4 CMEs 
in Group A occurred when the helicity increasing rate becomes much smaller (almost zero) after 
the monotonically increasing phase, and they are marked with the superscript 't ' in Table [TJ 
On the other hands, the time variations of AH for the 5 active regions in Group B are shown 
in Figure The source regions producing 17 CMEs in Group B presented a sign change of 
H before the CMEs: i.e., a noticeable increase in helicity of one sign was present for a while 
(6-45 hours) and then the CMEs occurred while relatively more helicity with the opposite sign 
was being injected in their source regions for the next few days. The main difference between 
Group A and Group B in the temporal variation of helicity injection is that the helicity sign 
reversal phase appears in Group B, but not in Group A. Note that we made the criterion for the 
helicity sign reversal shown in Group B as AH changes greater than 10 42 Mx 2 before and after 
the sign reversal phase starts. In the case that a sign change of H appeared but the criterion 
was not satisfied, we considered the variation of AH as fluctuation during the monotonically 
increasing phase of helicity injection shown in the active regions NOAA 8097, 8100, and 8693 
of Group A. 

We have also examined whether or not there is a specific trend in the temporal variation 
of $ in the CME source active regions during the period of the two characteristic helicity 
patterns shown in Group A and Group B. For doing that, from the time profiles of $ shown in 
Figures [2J [31 E2 and we calculated <& avg that is the average change rate of the total unsigned 
magnetic flux during the time period At between to and the CME occurrence time t\ in the 
active region under investigation: 

Mh) - $(i )] 



$ 



a eg 



At 



(?) 



In the case of the 30 CMEs in Group A, 19 CMEs (63 %) were preceded by an increase in $; 
i.e., <& avg > 0. However, the other 11 CMEs (37%) were preceded by a decrease in $; i.e., $ 



a eg 



-7 



< 0. For the 17 CMEs in Group B, 13 CMEs (76 %) occurred after an increase in $, while 4 
CMEs (24%) after a decrease in $. Refer to Tables __ and _0 for the detailed values of $(ti), 
^avgi and At of each event in Group A and Group B. We found that an increasing trend of $ 
before the occurrence of the CMEs is more common than a decreasing trend and this is more 
apparently shown in Group B than Group A. 

Now it would be interesting to know whether CME characteristics are different according 
to the two different patterns of helicity injection in the CME source regions. We therefore 
investigated a difference in the speed, acceleration, and occurrence trend of the CMEs between 
Group A and Group B. The CME speed and acceleration were adopted from the LASCO CME 



catalog (http://cdaw.gsfc.nasa.gov/CME_list). The linear speed in the catalog was taken 



as the representative CME speed in this study, which was determined from linear fit (also known 
as constant speed fit) to the height-time measurement of each of the CMEs within the LASCO 
C2 and C3 field of view (1.5-30 solar radii from the solar surface). The CME acceleration was 
also estimated from the height-time measurement using second-order polynomial fit. Note that 
there is a large uncertainty in some cases of the CME acceleration measurement with just three 
points in the height-time plot and they are marked with the superscript V in Tabled) 

We found that the CMEs in Group A have the speed range of 250-2460 kms -1 with the 
average (median) speed of 870(700) kms" 1 while those in Group B have the speed range of 500- 
2860 kms -1 with the average (median) speed of 1330(1150) kms -1 . The average CME speed of 
Group B is about 450 kms -1 faster than that of Group A. There is also a significant difference 
in the CME acceleration between Group A and Group B: the average acceleration is -24.4 ms -2 
and -6.3 ms -2 , respectively for Group A and Group B. These facts indicate that the CMEs in 
Group A have a relatively slow speed and rapid deceleration compared to those in Group B. 
Furthermore, the CMEs in Group A mainly consist of single events, while those in Group B 
tend to be successive events; i.e., 18 out of the total 23 active regions in Group A produced 
only one CME event during the helicity measurement period, while each of the total 5 active 
regions in Group B generated at least two CMEs. We examined soft X-ray characteristics of 
the CME-related flares in Group A and Group B using the 1-8 A GOES data. Both groups 
have a similar average flare duration, but they show a differe nce in total integrat ed flux F and 



flare impulsiveness / (i.e., peak flux/flare rise time; refer to lPearson et al.lll989l ): the average 
F (10 -5 Wm -2 min -1 ) is 1.5 and 2.3 and the average / (10 -5 Wm -2 min -1 ) is 0.7 and 1.1, 
respectively for Group A and Group B. This indicates that the CME-related flares in Group B 
are relatively more energetic and impulsive than those in Group A. In Table |2j the differences 
between Group A and Group B are summarized. 



-8 



3.2. Correlation bewteen CME Speed and Helicity Injection Rate 

From the time profile of H, a helicity parameter is defined to investigate its relationship 
with CME speed. We use the average helicity injection rate H avg which indicates the average 
amount of injected helicity per unit time into S during the time period At between t and t\. 

V' 1 Hit) 

tiavg - JTf , (» J 

where N is the total number of data points during At. Note that as a reference parameter, we 
consider $(ii), i.e., the total unsigned magnetic flux at t^. 



$&) = / \B n {x,t x )\dS. (9) 



s 



The two parameters H avg and 3>(ti) were calculated for the source regions producing the 
47 CMEs under investigation, and they were compared with the CME speed v as shown in 
Figure El The red and blue solid lines indicate the least-squares linear fits to the data points 
of Group A (red triangles) and Group B (blue diamonds) in each panel, respectively. The 
slope, intercept, and CC of the linear fits are also given in each panel. For Group A, we 
found that v has a very strong correlation (CC=0.85) with H avg , as well as a good correlation 
(CC=0.79) with $(ti). In case of Group B, v has a moderate correlation (CC=0.63) with H avg , 
while its correlation with $(£i) is weak (CC=0.45). This is because there is a high correlation 
(CC=0.83) between H avg and $(ti) for Group A (see Figure [7^) as previously reported that 
there is a good correlation be t ween helicity inje ction and total unsigned magnetic flux in active 
regions (jJeong fc Chad 120071 ; IPark et al.ll2010l ). but a poor correlation (CC=0.23) for Group 



B (as shown in Figure [7b). Note that the correlations between v and H avg for both Group 
A and Group B are statistically highly significant with > 99 % confidence level (from a two- 
tailed Student's t-test). In addition, we found that the linear fits to the data points of v 
(kms -1 ) vs. H avg (10 40 Mx 2 hr _1 ) are significantly different between Group A and Group B: 
i.e., v=14H avg +440 for Group A and v=128i/ at , 3 -|-700 for Group B. The slope of the linear 
fitted line for Group B is about nine times greater than that for Group A, and the intercept of 
the fitted line for Group B is ~250 km/s greater than that for Group A. 

We now have a question why there are considerable correlations between the helicity param- 
eter and the CME speed with a tendency that the lager injection of helicity a CME-productive 
active region achieves, the faster CME it produces. One possible answer for the question is as 
follows: the speed of CMEs depends on how much energy CMEs have, and the CME kinetic 
energy is supposed to be originated from free magnetic energy which is the energy deviation of 
the coronal magnetic field from its potential state. Assuming linear force-free magnetic fields, 
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Georgoulis fc LaBontd (120071 ) showed that there is linear dependence between free m agnetic 



energy and relative magn etic helicity from t he to tal magnetic energy formula given by iBerger 



f 19881 ). In the study of iRegnier fc Priestl (J2007I ). it was also found that there is a correla- 



tion between free magnetic energy and relative magnetic helicity estimated from reconstructed 
nonlinear force-free magnetic fields in four active regions even though the number of the active 
region samples is too small. Therefore, it might be reasonable to consider that helicity injection 
through the photosphere in a CME-productive active region is closely related to free magnetic 
energy build-up in the active- region corona for CME kinetic energy. 



4. Discussion and Conclusion 

47 CME events are examined to find a characteristic evolution pattern of helicity injection 
during a few days in their 28 source active regions. The main findings in this study are as 
follows: (1) there is always a significant helicity injection of 10 42 -10 44 Mx 2 through the active- 
region photosphere over a few days before the CMEs; (2) the CMEs under investigation are 
categorized into two different groups by two different helicity patterns which appeared for 
~0.5-4.5 before their occurrence. 30 CMEs in Group A occurred in 23 active regions in which 
monotonically increasing pattern with one sign of helicity injection is presented, whereas 5 
active regions producing 17 CMEs in Group B show a pattern of significant helicity injection 
followed by its sign reversal; and (3) the correlation between CME speed v and average helicity 
injection rate H avg in the CME-productive active regions is significant (CC=0.85) for Group A 
and moderate (CC=0.63) for Group B. 

In addition, the two CME groups classified by the two characteristic evolution patterns 
of helicity injection show different characteristics in kinematics, occurrence rate in a single 
active region, CME-related flare properties, and magnetic properties of their source regions 
as follows: (1) the average CME speed of Group B (1330 km s -1 ) is much faster than that of 
Group A (870 km s -1 ) and the CMEs in Group B (-6.3 ms -2 ) have a relatively slow deceleration 
compared to those in Group A (-24.4 ms -2 ); (2) a linear fit of v vs. H avg is quite different 
between Group A and Group B; (3) the CMEs in Group A tend to be single events, while 
those in Group B are inclined to be successive events; (4) soft X-ray flares related to the 
CMEs in Group B are relatively more energetic and impulsive than those in Group A; and (5) 
an increasing trend of magnetic flux in the CME source regions is more commonly shown in 
Group B than Group A. 

These differences may suggest that there are different pre-CME conditions for the two 
groups related to the two characteristic helicity patterns. We therefore try to understand how 
the two helicity patterns are involve with a precondition for CME initiation by comparing 
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them with temporal variations of helicity injection which are calculate d or expected from some 



numerical CME models. First, emerg i ng: tw i sted flux rope models (e.g..lTitov fc Demoulin 



Fan fc Gibson! 12004 ; iTorok fc Kliem! 120051 : Fan fe Gibson! 120071 : llsenberg fe Forbesl 120071 : [Fan 



1999 



2010l ) can be considered to explain the characteristic variation pattern o f A if shown in Group 
A (i.e., the monotonically increasing pattern with one sign of helicity). iFan &: Gibsonl (120041 ) 
showed in their numerical MHD simulation that as a twisted flux tube emerges gradually into 
a pre-existing coronal arcade during a pre-CME phase, the flux tube's helicity is continuously 
transported through the photosphere boundary into the corona. So we can verify that the 
time variation of photospheric helicity injection in the simulation is very similar to that in 
Group A (refe r to Fi gure 3c in Fan fc Gibsonl 120041 ) . In addition, the numerical simulations of 
Fan fc Gibsonl ( 120071 ) and iFanl (120101 ) indicated that (1) an emerging flux tube can settle into a 



phase of a steady quasi-static rise after the emergence of the flux tube is slowed down or stopped 
and (2) the quasi-static rising phase can be sustained for several hours until it erupts as a CME. 
We therefore conjecture that the phase of the nearly constant AH which appeared for several 
hours before the onset time of th e 4 CMEs in Group A may be related to the quasi-static stage . 
Second, some CME models (e.g-. lMikic et al.lll988[ Iwolfson fc Lowlll992l : lAntiochos et aDll999h 
which contain pre-CME dynamics associated with a steady shear (or twist) motion of magnetic 
field line footpoints on the photosphere can also explain the monotonically increasing pattern 
in Group A. In this case, magnetic flux emergence at the photosphere is not necessary to inject 
helicity, but helicity can be continuously injected via the steady photospheric shear (or twist) 
motion. 



On the other hand, we consider reconnection-fav ored emerging flux models (e.g.. lChen fc Shibata 



2000l : lArchontis et al.ll2007l ; lArchontis fc Hoodll2008l ) to explain the characteristic variation pat- 
tern of AH s hown in Group B (i.e., the helicity sign reversal pattern). The numerical MHD 
simulation of IChen fc Shibata] (120001 ) showed that a CME can be initiated by the localized re- 
connection between a pre-existing coronal field and a reconnection-favored emerging flux. And 
it is evident that this reconnection-favored flux can provide helicity injection of the opposite 
sign into an existing helicity system: hence, the phase of the opposite-signed helicity injection 
can be produced by a st eady and significant emergence of the reconnection-favored flux. In 
addition, as simulated by iKusano et al.l (j2003[ ). a relatively strong shear motion of underlying 
field footpoints in the reverse direction compared to that of overlying field footpoints can ex- 
plain the helicity sign reversal phase in Group B. Note also that there are some observational 
studies on spatial distributions of magnetic and electric current helicities in eruptive solar active 
regions, which may suggest a source region related to the helicity sign reversal. E.g., from the 
study of fractional el e ctric current helicity h c (i.e., B z ■ (V x B) z ) in 9 CME-associated active 



Wang et al.l (120041 ) reported that a key emerging flux region or a moving magnetic 



regions, 

feature near the main sunspot in the active regions brings up h c with a sign opposite to the 
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dominant sign of the main sunspot with a duration of a few days. iRomano et al.l (120111 ) studied 
a filament eruption in the active region NOAA 9682 and found that positive magnetic helicity 
was dominantly injected in the entire active region, while negative helicity was injected in local 
regions where the filament footpoints were located. 

In conclusion, we found that the two CME groups classified by the two characteristic 
helicity patterns have different pre-CME conditions and different CME characteristics. In 
addition, by comparing the observational helicity patterns with the expected patterns from 
some of numerical CME models, we presume that: the pre-CME condition of Group A is 
associated with an emergence of a twisted flux rope or a steady shear (or twist) motion of 
magnetic field line footpoints, while the pre-CME condition of Group B is involved with a 
reconnection-favored emerging flux or a reverse shear motion of magnetic field line footpoints. 
However, there are still some unanswered questions such as why there are significant differences 
between Group A and Group B in terms of CME kinematics (e.g., speed and acceleration) and 
why there is a fairly good correlation between H avg and $(ti) for the source regions in Group A, 
but a weak correlation for those in Group B. This helicity study should therefore be carefully 
checked with aspects shown in several other CME numerical simulations to further understand 
the physics underlying CME triggering mechanism and dynamics. And a future study will 
focus on the spatial distribution of helicity flux density and the evolution of detailed magnetic 
field structures in the investigated active regions to better understand CME initiation. 

Finally, we suggest that these characteristic variation patterns and helicity injection rate 
in CME-productive active regions can be used for the improvement of CME forecasting: (1) an 
early warning sign of CME occurrence could be given by the presence of a phase of monotonically 
increasing helicity as it is found that all the CMEs under investigation occur after significant 
helicity injection; (2) a warning sign for imminent CME occurrence could be also made when 
helicity injection rate becomes very slow or the opposite sign of helicity starts to be injected 
after the significant helicity injection in active regions; and (3) the speed of future CME event 
can be estimated by the statistical study of the correlation between the CME speed and the 
average helicity injection rate in the active regions. For doing this, more observational studies 
are therefore being carried out to check whether the two characteristic helicity pattern are 
shown in other CME-productive active regions and to improve the correlation between the 
CME speed and the helcity injection rate. 

The authors thank the SOHO/MDl team for the 96-minute full-disk photospheric mag- 
netogram data and the SOHO/LASCO team for the CME catalog which is generated and 
maintained at the CDAW Data Center by NASA and The Catholic University of America in 
cooperation with the Naval Research Laboratory. SOHO is a project of international coop- 
eration between ESA and NASA. We have made use of NASA's Astrophysics Data System 
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Table 1. List of 47 Selected CME Events 







CME 






Related Flare 


Source Region 




No. 


/ a 


PA/AW b 


v c 


a d 


e 


GOES 


NOAA 


Sunspot 


Group 




(deg) 


(kms -1 ) 


(ms - 2 ) 


Class 


No. 


Classification 




1 


1997-10-21 18:03 


Halo/360 


523 


-2.9 


17:00 


C3.3 


8097 


P 


A 


2 


1997-11-03 11:11 


232/122 


352 


-1.5 


09:03 


M1.4 


8100 


PlS 


At 


3 


1998-05-01 23:40 


Halo/360 


585 


8.0 


21:40 


C2.6 


8210 


7(5 


A 


1 


1998-05-02 05:31 


Halo/360 


542 


-1.4 


04:48 


C5.4 


8210 


/3 7 <5 


A 


5 


1998-05-02 14:06 


Halo/360 


938 


-28.8 


13:31 


XI. 1 


8210 


/3 7 <5 


A 


6 


1999-09-01 02:30 


188/283 


253 


-1.2* 


00:27 


C2.7 


8677 


P 


A 


7 


1999-09-13 09:30 


0/182 


898 


-23 


08:05 


C4.9 


8699 


P 


A 


8 


1999-09-13 17:31 


109/184 


444 


-8.7* 


16:30 


C2.6 


8693 


P 


A 


9 


1999-11-26 17:30 


228/145 


109 


6.0 


17:40 


C2.3 


8778 


P 


A 


10 


2000-01-18 17:54 


Halo/360 


739 


-7.1 


17:07 


M3.9 


8831 


P 


A 


11 


2000-01-28 20:12 


70/20 


429 


-2.8 


19:45 


C4.7 


8841 


P 


A 


12 


2000-05-10 20:06 


83/205 


611 


-15.5 


19:26 


C8.7 


8990 


P 


A 


13 


2000-06-06 15:54 


Halo/360 


1119 


1.5 


14:58 


X2.3 


9026 


PiS 


13 


11 


2000-06-07 16:30 


Halo/360 


812 


59.8* 


15:34 


XI. 2 


9026 


P-yS 


13 


15 


2000-07-14 10:54 


Halo/360 


1674 


-96.1* 


10:03 


X5.7 


9077 


PiS 


A 


16 


2000-07-25 03:30 


Halo/360 


528 


-5.8 


02:43 


M8.0 


9097 


Pi 


A 


17 


2000-08-09 16:30 


Halo/360 


702 


2.8 


15:19 


C2.3 


9114 


Pi 


A 


18 


2000-09-15 12:06 


249/235 


633 


-64.0* 


10:51 


C9.5 


9165 


ps 


A 


19 


2000-09-15 15:26 


217/210 


181 


-10.4* 


14:29 


M2.0 


9165 


ps 


A 


20 


2000-09-16 05:18 


Halo/360 


1215 


-12.3 


04:06 


M5.9 


9165 


Pi 


A 


21 


2000-10-02 03:50 


Halo/360 


525 


-4.9 


02:48 


C4.1 


9176 


Pi 


A 


22 


2000-10-09 23:50 


Halo/360 


798 


-9.8 


23:19 


C6.7 


9182 


P 


At 


23 


2000-11-23 23:54 


336/157 


690 


1.8 


23:18 


M1.0 


9236 


P 


13 


21 


2000-11-24 05:30 


Halo/360 


1289 


2.1 


04:55 


X2.0 


9236 


Pi 


13 


25 


2000-11-24 15:30 


Halo/360 


1245 


-3.3 


14:51 


X2.3 


9236 


Pi 


B 


26 


2000-11-24 22:06 


Halo/360 


1005 


-0.8 


21:43 


X1.8 


9236 


Pi 


13 


27 


2000-11-25 09:30 


Halo/360 


675 


-4.7 


09:06 


M3.5 


9236 


Pi 


13 


28 


2000-11-25 19:31 


Halo/360 


671 


-10.8 


18:33 


X1.9 


9236 


Pi 


13 


29 


2000-11-26 03:30 


259/118 


195 


-22.9* 


02:47 


M2.2 


9236 


Pi 


13 


30 


2000-11-26 17:06 


Halo/360 


980 


5.8 


16:34 


X4.0 


9236 


Pi 


13 


31 


2001-04-09 15:54 


Halo/360 


1192 


1.3 


15:20 


M7.9 


9415 


pyS 


13 


32 


2001-04-10 05:30 


Halo/360 


2411 


211.6* 


05:06 


X2.3 


9415 


0yS 


B 


33 


2001-09-17 08:54 


198/166 


1009 


-14.5 


08:18 


Ml. 5 


9616 


Pi 


At 


31 


2001-10-19 01:27 


Halo/360 


558 


-25.6 


00:47 


X1.6 


9661 


PiS 


A 


35 


2001-10-19 16:50 


Halo/360 


901 


-0.7 


16:13 


X1.6 


9661 


fSjS 


A 


36 


2002-07-15 20:30 


Halo/360 


1151 


-25.6 


19:59 


X3.0 


10030 


PlS 


13 


37 


2002-07-16 16:02 


Halo/360 


1636 


-41.0* 


- 


- 


10030 


pyS 


B 


38 


2002-08-16 12:30 


Halo/360 


1585 


-67.1 


11:32 


M5.2 


10069 


pyS 


A 


39 


2003-10-28 11:30 


Halo/360 


2459 


-105.2* 


09:51 


X17.2 


10486 


PlS 


A 


40 


2003-10-29 20:54 


Halo/360 


2029 


-146.5* 


20:37 


X10 


10486 


PlS 


A 


11 


2004-07-25 14:54 


Halo/360 


1333 


7.0 


14:19 


Ml.l 


10652 


PlS 


A 


12 


2004-11-06 01:31 


Halo/360 


818 


-81.5 


00:44 


M5.9 


10696 


PlS 


A 


13 


2004-11-07 16:54 


Halo/360 


1759 


-19.7 


15:42 


X2.0 


10696 


PlS 


At 


11 


2005-01-15 06:30 


Halo/360 


2049 


-30.7* 


05:54 


M8.6 


10720 


ps 


B 



-16- 



Table 1 — Continued 







CME 






Related Flare 


Source Region 




No. 


f„ a 


PA/AW b 


v c 


a d 


GOES 
Class 


NOAA 


Sunspot 


Group 




(dcg) 


(kms -1 ) 


(ms - 2 ) 


No. 


Classification 




45 


2005-01-15 23:06 


Halo/360 


2861 


-127.4* 


22:25 X2.6 


10720 


/3<5 


13 


46 


2005-01-17 09:54 


Halo/360 


2547 


-159.1* 


06:59 X3.8 


10720 


06 


13 


47 


2006-07-04 21:30 


199/102 


308 


1.6 


19:06 C1.4 


10898 


P 


A 



Note. — 

a First appearance time in the LASCO/C2 FOV. 

b Position Angle/Angular Width given by the LASCO CME catalog. 

c Linear speed given by the LASCO CME catalog. 

d Acceleration given by the LASCO CME catalog. The superscript '*' indicates that acceleration is uncertain due 
to either poor height measurement or a small number of height-time measurements. 

° GOES soft X-ray flare start time. 

The CMEs under investigation are classified into two groups by two patterns of helicity injection in their source 
regions: a pattern of monotonically increasing helicity in Group A and a pattern of sign reversal of helicity in Group 
B. The superscript 'f indicates the CME events in Group A which occur during a phase of relatively constant helicity 
after a significant amount of helicity injection. For the detailed description of Groups A and B, refer to § 13.11 
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Table 2. Comparison of Different Characteristics between Group A and Group B 



Group 



CME 
Number 



b CME CME-related Flare 

lavg Time. J? c T d 

lype ^avg Javg 



Flux Pattern 
®avg > $avg < 



Helicity Pattern 



A 30 870 -24.4 Single 1.5 0.7 63% 37% Monotonic increase 

B 17 1330 -6.3 Successive 2.3 1.1 76% 24% Sign reversal 

a Average velocity (kms -1 ). 

b Average acceleration (ms~ 2 ). 

c Average total time-integrated flux measured in the 1 to 8 A by the GOES satellite (10 — 1 Jm -2 ). 

d Average flare impulsiveness which is determined from the ratio of peak flux to flare rise time (10 — 5 Wm~ 2 mill -1 ). 

c <&avg is average change rate of total unsigned magnetic flux in a CME source region. 
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Table 3. Helicity Injection in 23 Active Regions Producing 30 CMEs in Group A 





CME 




Source Re 


gion 






No. a 


V 


a 


flavg 


*(ti) c 




<i> d 


Ar c 




(kms- 1 ) 


(ms~ 2 ) 


(lO^Mx 2 !!!-- 1 ) 


(10 20 Mx) 


(10 20 


Mxday -1 ) 


(day) 


1 


523 


-2.9 


0.1 


180 




-20 


2.1 


2t 


352 


-1.5 


5.1 


530 




61 


3.3 


3 


585 


8.0 


11.5 


460 




14 


3.0 


1 


512 


-1.4 


12.9 


460 




15 


3.2 


5 


038 


-28.8 


15.9 


480 




19 


3.6 


6 


253 


-1.2* 


0.7 


120 




-18 


3.1 


7 


898 


-23 


8.1 


300 




83 


1.7 


8 


444 


-8.7* 


0.2 


100 




-13 


2.5 


9 


409 


6.0 


8.8 


390 




30 


2.1 


10 


73,0 


-7.1 


23.5 


560 




-19 


1.7 


11 


120 


-2.8 


3.6 


210 




7 


2.8 


12 


641 


-15.5 


9.7 


360 




97 


0.5 


15 


1674 


-96.1* 


59.4 


710 




-15 


2.5 


16 


528 


-5.8 


14.8 


440 




3 


2.7 


17 


702 


2.8 


25.9 


400 




23 


2.7 


18 


633 


-64.0* 


23.8 


450 




81 


2.1 


19 


181 


-10.4* 


25.8 


150 




78 


1.7 


20 


1215 


-12.3 


28.5 


460 




66 


2.8 


21 


525 


-4.9 


31 


530 




-9 


1.0 


09t 


798 


-9.8 


2.1 


320 




7 


3.5 


33+ 


1009 


-14.5 


8.4 


620 




-11 


1.2 


31 


558 


-25.6 


30.2 


810 




9 


2.1 


35 


901 


-0.7 


40.6 


810 




7 


2.8 


38 


1585 


-67.1 


57.8 


1010 




80 


1.1 


30 


2459 


-105.2* 


128.5 


1440 




-12 


1.0 


10 


2029 


-146.5* 


104.8 


1470 




-3 


2.1 


41 


1333 


7.0 


70.6 


1110 




-10 


4.6 


42 


818 


-81.5 


90.7 


450 




77 


1.1 


43+ 


1759 


-19.7 


68.5 


5 10 




61 


2.7 


47 


308 


1.6 


9.9 


470 




-1 


2.5 



a The ID numbers here are corresponding to those in Table [T] 

b Average of helicity injection rate in the active region under investigation from the start 
time of helicity measurement to to the CME occurrence time ti . 

c Total unsigned magnetic flux in the active region under investigation at the CME occur- 
rence time t\. 

d The average change rate of the total unsigned magnetic flux during the time period At 

e The time period between the start time of helicity measurement and the CME occurrence 
time. 
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Table 4. Helicity Injection in 5 Active Regions Producing 17 CMEs in Group B 



CME Source Region 

No. a V a Havg *(tl) &avg At 

(kms- 1 ) (ms" 2 ) (10 40 Mx 2 hr- 1 ) (10 20 Mx) (10 20 Mxday" 1 ) (day) 



13 


1119 


1.5 


2.4 


830 


-51 


1.1 


14 


842 


58.8 


0.9 


780 


-51 


2.1 


23 


690 


1.8 


7.9 


600 


22 


1.7 


21 


1289 


2.1 


5.8 


610 


2.5 


1.9 


25 


1245 


-3.3 


2.2 


610 


35 


2.3 


26 


1005 


-0.8 


0.5 


660 


38 


2.6 


27 


675 


-4.7 


1.0 


680 


10 


3.1 


28 


671 


-10.8 


4.7 


710 


11 


3.5 


29 


195 


-22.9 


2.5 


710 


10 


3.8 


30 


980 


5.8 


0.5 


7.50 


12 


1.1 


31 


1192 


1.3 


8.3 


660 


-8 


2.5 


.32 


2411 


211.6 


6.1 


660 


-9 


3.0 


36 


1151 


-25.6 


7.2 


900 


69 


2.3 


3V 


1636 


-41.0* 


3.5 


980 


75 


3.1 


11 


2049 


-30.7* 


5.0 


810 


103 


1.7 


15 


2861 


-127.4 


11.1 


850 


90 


2.1 


46 


2547 


-159.1 


12.0 


940 


79 


3.8 





The ID numbers here are corresponding to those in Table [T] 




Fig. 1.— Solar active region NOAA 10720 at 14:23 UT on 2005 January 16. Panels (a) and 
(b) show ^lct (red/yellow arrows on positive/negative B n ) and A p (red/yellow arrows on 
positive/negative B n ) superposed on the grayscale image of B n derived from the MDI line-of- 
sight magnetogram, respectively. In Panel (c), Ga map is presented in grayscale. Note that 
the saturation level of \Ga\ is set as lxlO 40 Mx 2 hr -1 Mm -2 for purpose of display visibility. 
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Fig. 2. — Time variation of helicity injection AH (black crosses) and total unsigned magnetic 
flux $ (blue diamonds) for 8 active regions in Group A. The active regions in Group A show a 
monotonically increasing pattern of helicity for a few days. In each panel, the vertical red lines 
indicate the times when the CMEs originating from the 8 active regions first appeared in the 
LASCO/C2 FOV. 
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Fig. 3. — Same as in Figure [2} but for additional 8 active regions in Group A. 
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Fig. 4. — Same as in Figure [2} but for additional 7 active regions in Group A. 
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Fig. 5. — Same as in Figure [2j but for 5 active regions in Group B indicating a pattern of 
significant helicity injection followed by its sign reversal. A total of 17 CMEs occurred from the 
5 active regions during the period when the helicity injection rate in the active regions started 
to reverse its sign. 
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Fig. 6. — The speed v of 47 CMEs (30 CMEs in Group A marked by red triangles and 17 
CMEs in Group B marked by blue diamonds) is plotted against the two magnetic parameters 
of (a) the average helicity injection rate H avg and (b) the total unsigned magnetic flux $(ii). In 
each panel, the red and blue solid lines indicate the least-square linear fits to the data points of 
Group A and Group B, respectively, and the slope and intercept of the fitted lines are specified 
with the linear correlation coefficient (CC). The correlation between v and H avg is very strong 
(CC=0.85) for Group A and moderate (CC=0.63) for Group B. See Table [3] for the detailed 
information of the 47 CMEs. 
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Fig. 7. — The total unsigned flux $(ti) vs. the average helicity injection rate H avg for (a) 30 
CMEs in Group A and (b) 17 CMEs in Group B. The least-square linear fit (solid line) and 
the linear correlation coefficient (CC) are presented in each panel. A very high correlation 
between H avg and $(ii) is found for Group A (CC=0.83), but a poor correlation for Group B 
(CC=0.23). 



